CLNS 07/1998 
CLEO 07-07 

Measurement of the total hadronic cross section in e + e 
annihilation below 10.56 GeV 

D. Besson, 1 T. K. Pedlar, 2 D. Cronin-Hennessy, 3 K. Y. Gao, 3 J. Hietala, 3 T. Klein, 3 
Y. Kubota, 3 B. W. Lang, 3 R. Poling, 3 A. W. Scott, 3 A. Smith, 3 P. Zweber, 3 S. Dobbs, 4 
Z. Metreveli, 4 K. K. Seth, 4 A. Tomaradze, 4 J. Ernst, 5 K. M. Ecklund, 6 H. Severini, 7 

S. A. Dytman, 8 W. Love, 8 V. Savinov, 8 O. Aquines, 9 A. Lopez, 9 S. Mehrabyan, 9 
H. Mendez, 9 J. Ramirez, 9 G. S. Huang, 10 D. H. Miller, 10 V. Pavlunin, 10 B. Sanghi, 10 
; I. P. J. Shipsey, 10 B. Xin, 10 G. S. Adams, 11 M. Anderson, 11 J. P. Cummings, 11 I. Danko, 11 

D. Hu, 11 B. Moziak, 11 J. Napolitano, 11 Q. He, 12 J. Insler, 12 H. Muramatsu, 12 C. S. Park, 12 
^ E. H. Thorndike, 12 F. Yang, 12 M. Artuso, 13 S. Blusk, 13 J. Butt, 13 J. Li, 13 N. Menaa, 13 

R. Mountain, 13 S. Nisar, 13 K. Randrianarivony, 13 R. Sia, 13 T. Skwarnicki, 13 S. Stone, 13 
J. C. Wang, 13 K. Zhang, 13 G. Bonvicini, 14 D. Cinabro, 14 M. Dubrovin, 14 A. Lincoln, 14 
D. M. Asner, 15 K. W. Edwards, 15 P. Naik, 15 R. A. Briere, 16 T. Ferguson, 16 G. Tatishvili, 16 
^ I H. Vogel, 16 M. E. Watkins, 16 J. L. Rosner, 17 N. E. Adam, 18 J. P. Alexander, 18 

K. Berkelman, 18 D. G. Cassel, 18 J. E. Duboscq, 18 R. Ehrlich, 18 L. Fields, 18 R. S. Galik, 18 

L. Gibbons, 18 R. Gray, 18 S. W. Gray, 18 D. L. Hartill, 18 B. K. Heltsley, 18 D. Hertz, 18 
C. D. Jones, 18 J. Kandaswamy, 18 D. L. Kreinick, 18 V. E. Kuznetsov, 18 H. Mahlke-Kriiger, 18 
D. Mohapatra, 18 P. U. E. Onyisi, 18 J. R. Patterson, 18 D. Peterson, 18 J. Pivarski, 18 
D. Riley, 18 A. Ryd, 18 A. J. Sadoff, 18 H. Schwarthoff, 18 X. Shi, 18 S. Stroiney, 18 
W. M. Sun, 18 T. Wilksen, 18 S. B. Athar, 19 R. Patel, 19 V. Potlia, 19 J. Yelton, 19 
P. Rubin, 20 C. Cawlfield, 21 B. I. Eisenstein, 21 I. Karliner, 21 D. Kim, 21 N. Lowrey, 21 
m i M. Selen, 21 E. J. White, 21 J. Wiss, 21 R. E. Mitchell, 22 and M. R. Shepherd 



00 
(N 



O 



22 



(CLEO Collaboration) 

\q . 1 University of Kansas, Lawrence, Kansas 66045 

2 Luther College, Decorah, Iowa 52101 
3 University of Minnesota, Minneapolis, Minnesota 55455 
^ Northwestern University, Evanston, Illinois 60208 
5 State University of New York at Albany, Albany, New York 12222 
6 State University of New York at Buffalo, Buffalo, New York 14260 
7 University of Oklahoma, Norman, Oklahoma 73019 

8 University of Pittsburgh, Pittsburgh, Pennsylvania 15260 

9 University of Puerto Rico, Mayaguez, Puerto Rico 00681 
10 Purdue University, West Lafayette, Indiana 47907 

11 Rensselaer Polytechnic Institute, Troy, New York 12180 
12 University of Rochester, Rochester, New York 14627 
13 Syracuse University, Syracuse, New York 13244 
* Wayne State University, Detroit, Michigan 48202 
15 Carleton University, Ottawa, Ontario, Canada K1S 5B6 
16 Carnegie Mellon University, Pittsburgh, Pennsylvania 15213 
17 Enrico Fermi Institute, University of Chicago, Chicago, Illinois 60637 
18 Cornell University, Ithaca, New York 14853 

19 University of Florida, Gainesville, Florida 32611 

20 George Mason University, Fairfax, Virginia 22030 
21 University of Illinois, Urbana- Champaign, Illinois 61801 



1 



22 Indiana University, Bloomington, Indiana 4^405 
Abstract 

Using the CLEO III detector, we measure absolute cross sections for e + e~ — > hadrons at seven 
center-of-mass energies between 6.964 and 10.538 GeV. The values of R, the ratio of hadronic and 
muon pair production cross sections, are determined within 2% total r.m.s. uncertainty. 
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I. INTRODUCTION 



The total cross section for hadron production in electron-positron annihilation is often 
expressed in terms of the ratio R(s) of the radiation-corrected measured hadronic cross 
section to the calculated lowest-order cross section for muon pair production. That is, 



R(s) = <7 (e + e — > hadrons) /cr (e + e — > 



(1) 



where cr (e + e~ — > = 4:ira 2 /3s . In the simplest form of the quark model R(s) is Rq(s), 

the sum of the squares of the charges for the quark species allowed by energy conservation, 
including a factor of iV c = 3 for the colors. For energies between cc and bb thresholds, 
Rq = 3S udsc Qj = 10/3, where Qf represents the charge of the quark f — u,d, s, c. 

The strong interactions of the quarks can be accounted for in a perturbation theory 
expansion in the QCD coupling az s (s). In the modified minimal-subtraction scheme at the 
four- loop level [l[ : 



R(s) = R 



71 



71 



+ C3 ln) +0 ^ S )) 



(2) 



For the four-flavor case the coefficients are d = 1, C 2 = 1.525, and C 3 = -11.686 Q, Q. 
Asymptotic freedom in QCD implies that the coupling a s is not a constant but a smoothly 
decreasing function of s. The aim of the present experiment is to use measurements of R(s) 
to make an accurate test of this relation in the 48.5 to 111.0 GeV 2 region of s, where u, d, s, 
and c quarks are produced. 

We use data collected at seven center-of-mass energies at the Cornell Electron Storage 
Ring (CESR). Electrons and positrons are stored at equal energies in oppositely rotating 
orbits in the Cornell Electron Storage Ring (CESR), and they collide at a 2 mrad angle at 
the center of the CLEO III detector [3j . Charged tracks are reconstructed in a 4-layer silicon 
strip detector and a 47-layer wire drift chamber, and their momenta are determined from 
their radii of curvature in a 1.5 T magnetic field created by a superconducting solenoid with 
its axis aligned along the average of the two colliding beam lines. The Csl scintillator shower 
calorimeter forms a cylindrical barrel around the tracking volume, covering polar angles 9 
with respect to the beam axis with | cos 9 |< 0.85, supplemented by endcaps extending the 
range to | cos 9 |< 0.93. Electromagnetic showers are detected with an r.m.s. resolution of 
about 75 MeV at 5 GeV. 



II. ANALYSIS OVERVIEW 

To obtain the continuum hadronic annihilation cross section referred to in Eqs. ([HE]), we 
must correct the observed yield for detection efficiency, for initial state radiation including 
its influence on efficiency, for other higher order QED effects, and for contributions from 
the radiative tails of lower energy resonances. We write the observed continuum hadronic 
cross section as the sum of a sv from soft photon radiation and virtual higher order diagrams, 
plus o"hard from hard photon radiation, plus <r res from radiative tails of lower mass qq bound 
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states QM, B0 : 



O"obs(s) = OsvO) + (T hard (s) + CT res (s). (3) 

The sum of the soft and virtual term can be expressed as 

a sv (s) = e(0)a (s)5 S v, (4) 

where e(0) is the detection efficiency for events with no initial state radiation, <jq is the 
desired corrected hadronic Born cross section, and S sv = 5 + 5 vp accounts for soft photon 
emission and hadronic and leptonic vacuum polarization: 



2a / 3 , s 7r 2 > 
5 = — -In— + — -1 , 5 
7r \ 4 mi 6 ' 



and 



where / = e,/i,T,h. A more accurate approximation of hadronic vacuum polarization can 
be found in reference 

The hard photon piece involves an integral of the continuum cross section over all radiated 
energies 0: 

O r hard(s) = O- (s)J hard , (7) 

with 

where /c = (s — s')/s is the normalized radiated energy, -v/? is the center of mass energy 
after the initial state radiation, and the number of equivalent radiation lengths is 

*=^fln-L-l). (9) 

Ti \ m 2 e J 

Before making corrections to the observed continuum hadronic cross section, we first 
have to subtract the contributions from lower energy qq resonances excited through the 
radiative process e + e~ — > 'jqq (the last term in Eq. Q). For each contributing resonance 
the corresponding term in cx res is given by the product of the detection efficiency e(k), the 
non-radiative cross section energy integral (given in terms of the resonance coupling T ee to 
e + e~, the resonance hadronic branching ratio -Bhad, and the resonance mass M), and the 
radiative kernel, F(k,s): 



Cres(s) 



T,<k)(^r ee B h ^F(k,s), (10) 
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where the sum runs over all known low energy resonances which contributes observed 
hadronic cross section. In lowest order the kernel can be written as 



with k = (s — s')/s = 1 — M 2 / s. A higher order representation of F(k, s) can be found in 



reference 10 



The goal of the measurement, the Born cross section, is derived from Eqs. (131), (HI), and 



/ x cTobsls) - a res (s) 
e(0)d sv + i hard 

It is obtained from the observed hadronic cross section, cr b s = (-/Vhad — -^Vbkg)/(f trigger x / Cdt), 
where iVhad is the number of detected events passing hadronic event selection criteria, iVbkg 
is the number of background events passing the same requirements, e tr igger is the trigger 
efficiency correction, and / Cdt is the integrated luminosity. The background includes 
e, /x, and r pairs, virtual photon-photon (two-photon) interactions (e + e~ — > e + e~Yl* ~^ 
e + e~ hadrons), beam-gas (beam particle interactions with residual gas molecules) and beam- 
wall interactions (beam particle interacts with the beam pipe material), and cosmic rays. 



III. DATA SAMPLES 

The center-of-mass energies y/s and integrated luminosities / Cdt used here are listed in 
Tabled 

TABLE I: Center-of-mass energies and integrated luminosities. 



Energy 


Luminosity 


GeV 


pb" 1 


10.538 


904.50 ± 0.30 ± 9.00 


10.330 


149.80 ±0.10 ±1.60 


9.996 


432.60 ± 0.20 ± 4.80 


9.432 


183.00 ±0.10 ±2.00 


8.380 


6.78 ± 0.02 ± 0.06 


7.380 


8.48 ± 0.02 ± 0.07 


6.964 


2.52 ±0.01 ±0.02 



The luminosities for these data were calculated from measured rates of three different 
QED processes, e + e~ — > e + e~, 77, and yU + yU~, using predicted cross sections from the 



Babayaga MC generator [11|], with selection criteria similar to those in Ref. [12]. Particular 
care was taken to include the effects of T(nS) — > / + / _ , where I stands for lepton: e or /x, 
from either direct production or with an intermediate radiative return [10J. The resulting 
luminosities from all three reactions were found to be consistent within their systematic 
uncertainties of 1.8%, 1.2%, and 1.2% for the 77, e + e~, and channels, respectively, 
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and a single value at each energy was obtained by combining the three values according to 
their relative total uncertainties. The total relative uncertainty on the combined luminosity 
varied from 0.9-1.1%, depending upon energy. 

To study efficiencies and backgrounds we used GEANT-based Monte Carlo simulations 



13j, Jetset 7.3 [14( for hadronic events, and KORALB 15] for r pairs. 



IV. EVENT SELECTION AND DETECTION EFFICIENCIES 

The detector trigger used in the measurement of the hadronic event rate is the logical 
"or" of the following two sets of requirements: 1) at least three axial-projection charged 
tracks (tracks which are detected in the inner 16 layers of the drift chamber) and at least 
one low-threshold electromagnetic shower; 2) at least two stereo-projection charged tracks 
(tracks which are detected in the outer 31 layers of the drift chamber) and either at least two 
low-threshold showers or at least one medium-threshold shower. The low (medium) trigger 
threshold for showers is approximately 150 MeV (750 MeV). We evaluate the efficiency by 
comparing rates (using the hadronic selection criteria described in the next section) with 
those obtained in special runs with a much looser detection trigger, namely requiring two or 
more axial-projection charged tracks. Table [Til shows the resulting trigger efficiency values, 
^trigger, for the various energy points. The systematic uncertainties are conservatively taken 
as the deviations from 100% efficiency. 

TABLE II: Trigger efficiency correction cogger, and net hadronic event detection efficiency e(0) 
evaluated for the case of zero initial state radiation energy (k = 0). The uncertainties are statistical 
only. 



Energy 


^trigger 


6(0) 


GeV 


% 


% 


10.538 


99.91 ±0.01 


87.43 ± 0.01 


10.330 


99.91 ±0.01 


87.42 ± 0.03 


9.996 


99.89 ± 0.01 


87.13 ±0.03 


9.432 


99.92 ±0.01 


86.05 ± 0.04 


8.380 


99.88 ± 0.02 


84.71 ±0.10 


7.380 


99.87 ±0.02 


82.95 ±0.10 


6.964 


99.81 ± 0.03 


82.07 ±0.14 



In order to suppress backgrounds from events other than e + e~ — > hadrons, we apply 
selection requirements (cuts) to individual tracks and showers as well as to entire events. 
We now describe these cuts and explain how we estimate their efficiency and associated 
uncertainties. 

Table III 1 1 lists the requirements for accepting individual tracks and showers. We accept 
all tracks which have \ 2 P er degree of freedom (x 2 /NDF) less than 100. Consistency with 
the beam collision point is enforced by the cut on do, the distance of closest approach of the 
reconstructed track relative to the beam axis, and on zo, the distance between that point 
and the average collision point on the beam axis. The hit fraction is the ratio of the number 
of detected and expected hits when a track passes through the drift chamber, and the 9 is 
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TABLE III: Cuts used to select tracks and showers. 



Variable Allowed range 

X 2 /NDF < 100.0 

hit fraction (0.5, 1.2) 

| do I < 3.0 cm 

| Zq I < 18.0 cm 

error of zq < 25.0 cm 

| cot (6>) | < 3.0424 

error of cot (8) < 0.50 

-ftrack/^beam (0.01,1.5) 

^shower /-^beam 0.01 



the polar angle of the track relative to the beam axis. We require that each track candidate 
carry at least 1% of the beam momentum. To avoid poorly measured tracks, we require that 
the ratio of the track momentum to the beam momentum be less than 1.5. We require that 
each accepted shower carry at least 1% of the beam energy and that no track be associated 
with it. 

The event selection criteria are listed in Table IIVI and their description are provided in 
the following paragraphs. 



TABLE IV: Hadronic event selection criteria. 



Variable Allowed range 

| -^vertex | < 6.0 Cm 

-E'vis/(2-E'beam) > 0.5 

| Pr ss /E vis | < 0.3 

H 2 /H < 0.9 

£ ca i/(2£ beam ) (0.15,0.9) 

£™ ax /i?beam < 0.8 

AfchargedTrack > 4 



The event Z veItex is the weighted average of the z^s of the charged tracks, and weights 
are defined as l/of , where a zo is the uncertainty of the z . The purpose of this cut was 
to suppress background due to beam-gas, beam-wall collisions and from cosmic rays. The 
measured R value was sensitive at the level of only 0.16% to the cut over the range from 4 
to 10 cm. 

The visible energy E V - 1S is defined as the sum of charged particle energies computed from 
the track momenta assuming pion masses, plus the neutral shower energies measured in the 
electromagnetic calorimeter. We normalize the visible energy to the center-of-mass energy 
so that we can use the same cut for all runs. The visible energy cut is designed to suppress 
background from virtual photon-photon collisions and from beam-gas interactions. Figure [1] 
shows the normalized visible energy distribution at 9.996 GeV. We varied the cut in scaled 
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4180507-001 




FIG. 1: Visible energy distribution at 9.996 GeV, normalized to center-of-mass energy, when all 
cuts are applied except the one shown. In this as well as in all other figures, arrows indicate cut 
values. 

visible energy (E vis / (2£ , b ea m)) from its nominal level of 0.5 down to 0.3 and to 0.4 to estimate 
our sensitivity to the two-photon background. We also varied the cut upwards to 0.6 and 
to 0.7 to estimate the discrepancy between data and Monte Carlo. 

The background from two-photon and beam-gas events was also suppressed by the cut 
on the ratio | P™ ss /E vis | of the ^-component of the event missing momentum and visible 
energy (see Fig. [2]). The R results varied by less than 0.65% over the cut range 0.2 to 0.4. 

4180507-002 

10 5 
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d 1 ° 4 

CD 

w 103 



-0.8 -0.6 -0.4 -0.2 -0 0.2 0.4 0.6 0.8 
pmiss/p 

FIG. 2: The ratio of the z-component of the event missing momentum to event visible energy 
P™ iss /E vis at 10.538 GeV. 

The typical signature of a QED event (e, /i, or r pair) is low charged track multiplicity 
which would fail the four-track requirement, although higher multiplicity decays of taus are 
likely to pass. We rely on calorimeter information to suppress the remaining background 
further. The distribution of the ratio of calorimeter energy E ca ^, which includes isolated 




S 



showers as well as track matched showers, to the center-of-mass energy, is plotted in Fig. [3j It 
shows a contribution at low values from tau-pairs and, at E ca \/ {2E\ )Cara ) = 1, a contribution 
from Bhabha-scattered electron pairs. The cuts eliminated both extremes. Varying the 
upper cut from 0.8 to oo made less than 0.49% difference in R values. 
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FIG. 3: The ratio of calorimeter energy of charged and neutral particles to 2Ei ocam at 10.538 GeV. 



Figure H] shows the distribution of the ratio of Fox- Wolfram moments H^/Hq 16| . The 
QED process e + e~ — ► l + l~, where I — e, \i or r, results in a high value of this ratio. Cut 
variation above 0.8 produced a variation in R of less than 0.54%. 
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FIG. 4: Distribution of the ratio of Fox- Wolfram moments H2/H0 at 6.964 GeV. 



We also studied the ratio of the energy of the most energetic unmatched shower to the 
beam energy, shown in Fig. [5j To suppress production of hadronic events at low center-of- 
mass energies through initial state radiation, we cut at 0.80. Varying the cut from 0.6 to 
infinity changed the value of R less than 0.40%. 

The charged track multiplicity cut also suppressed QED events in our hadronic sample. 
The multiplicity distribution (see Fig. [6]) showed a significant sensitivity of the accepted 
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FIG. 5: The ratio of the most energetic shower energy to the beam energy at 8.380 GeV. 

rate to the chosen cut level. However, there is a discrepancy between data and the Monte 
Carlo simulation. The discrepancy at lower multiplicity values can be explained by the 
contribution of QED events; the discrepancy at higher multiplicity values can be explained 
by inaccuracies at the Monte Carlo generator level. 
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FIG. 6: Charged particle multiplicity distribution at 10.330 GeV. 



The selection criteria defined in Tables III 1 1 and [IV] were used to determine the number 
of hadronic events and the hadronic detection efficiencies for the various incident beam 
energies, using our data and Monte Carlo samples (see Tables HT1 and IVl) . We performed the 
efficiency calculation using weights assigned to each continuum multiplicity value instead 
of relying on the continuum Monte Carlo prediction for the multiplicity distribution. The 
weights were assigned by comparing data multiplicity with continuum and r-pair Monte 
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Carlo samples: 



w,. 



yydata _ yyr 
n n 



-pair 



(13) 



where w n is the weight which is assigned to a continuum event with multiplicity n; iV^ ata , 
iV r_pair , and N^ ont are the number of events in the data sample, r-pair, and continuum 
Monte Carlo samples with multiplicity n, respectively. Since low multiplicity values are 
contaminated by e + e~ and events, we assign a weight of 1.0 for multiplicities less than 

three. 

Figure [7] shows the fc-spectra, relative decrease in squared center-of-mass energy due 
to initial state radiation, The step near k = 0.84 represents the cc threshold. Table HT1 
summarizes the no-radiation hadronic event detection efficiencies e(0), and Fig. [8] shows 
how the hadronic event detection efficiency varies with k. 
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FIG. 7: k distribution of the Monte Carlo continuum sample at 10.538 GeV. 



V. BACKGROUND 

Sources of background were suppressed with the help of our hadronic event selection 
criteria. We deal with the remaining level of background events from two-photon, beam- 
gas, and beam-wall interactions, as well as from r-pair, /i-pair, and Bhabha production in 
the following. Figure M shows the net charge distribution of hadronic events that pass the 
hadronic event selection criteria. The Monte Carlo (continuum + r-pair) describes our data 
well, however, there are small discrepancies in the tails (Figs. 1-6). These discrepancies 
may be related to the beam-gas or beam-wall interactions not included in our Monte Carlo 
samples. We estimated and subtracted the number of r-pair events using the known QED 
cross section and the KORALB [l5| Monte Carlo simulation. 

The two-photon interactions, which are the dominant source of background after r-pairs, 
were studied with a dedicated Monte Carlo that generated e + e~ — > e + e~ + hadrons events. 
We found that the portion of two-photon events which survive our selection criteria with a 
loose cut E vis /2E heam > 0.2 (our nominal cut is 0.5) changes from 1.52% at 10.538 GeV to 
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FIG. 8: Hadronic event detection efficiency as a function of k at 9.996 GeV. 
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FIG. 9: Net charge distribution of events at 9.996 GeV. 

0.5% at 6.964 GeV. For this reason, instead of subtracting the two-photon background, we 
have varied the E vis /2E^ eam cut to assign an additional contribution to the event selection 
systematic uncertainty in R. 

VI. RESONANCE EFFECTS AND RADIATIVE CORRECTION 



The beams radiate energy before the annihilation interaction and can then create cc 
and bb bound state resonances. When the resonance decays to hadrons, our observed cross 
section increases. For the purposes of our R measurement, these contributions are sources 
of background and must be subtracted. The following resonances were considered: T(3S), 
T(2S), T(1S), '0(4160), -0(3770), 0(2S), and J/if). To take into account the resonance 
contributions we use Eqs. ( TTOT) and (j3J). For this purpose we used the CLEO-measured 
dielectron widths of T(1S), T(2S), and T(3S) from Ref. [3] and dielectron widths of J/ip 
and -0(2S) from Ref. [19J. We used dielectron widths of -0(3770) and -0(4160) from the 
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PDG 2(j. Table M shows the total cross section of contributing resonances and the ratio of 
the resonance contribution to the observed cross section for each energy point. We used a 
sample of low-background charged dipion transition events in data to check our calculated 
quarkonium contributions to our observed hadronic sample. Selecting T(3S), T(2S), and 
■0(2S) resonance production through initial state radiation from 10.538 GeV, we studied the 
following transitions: T(3S) -> 7r+7r-T(lS), T(2S) -> 7r+7r-T(lS), and ^(2S) -> tt+tt- J/i/>, 
where T(2S), T(1S), and J /if) decay to to estimate resonance contributions at our 

highest-luminosity point, 10.538 GeV. We found that our experimental results agreed with 
the theoretical expectation for these modes within 1.5 standard deviations. 



TABLE V: Number of hadronic events after background subtraction, radiation and efficiency cor- 
rections, total and relative contribution of resonances through radiative return, and interference 
effect. 



Energy 


Hadronic 


e(0)5 sv + 4ard 


Ores 


0-res/o-obs Interfer- 


GeV 


events 




nb 


on 


ence % 


10.538 


2426337 


0.935 


0.056 


2.08 




10.330 


403743 


0.935 


0.044 


1.62 


0.14 


9.996 


1988179 


0.931 


0.033 


1.14 


0.22 


9.432 


573516 


0.919 


0.001 


0.03 


0.51 


8.380 


26705 


0.901 


0.003 


0.09 




7.380 


42204 


0.879 


0.009 


0.14 




6.964 


14508 


0.867 


0.009 


0.17 





Since three runs were at energies just below the Upsilon resonance masses, interference 
between the continuum and the resonance can affect the measured continuum cross section. 
The magnitude and sign of this effect depends on the difference between the center-of-mass 
energy and the resonance mass. The CESR beam energy spread, which is of order 4 MeV, 
convolves the interference effects from different center-of-mass energies, yielding an overall 
negative effect on the measured cross section below resonance. In the last column of Table IVl 
are listed calculated corrections to the continuum cross section to compensate upward for 
the effect of interference and beam energy spread [l8| . 

We applied the procedure described in Section fll] by evaluating the e(0)<5 sv + /hard term of 
Eq. (j!2p to calculate the radiative correction for each incident energy point. The efficiencies 
at different values of k (see Fig. [8]) from continuum Monte Carlo have been used to calculate 
-fhard- Table |V] shows the radiation and efficiency factor [see Eq. (1121) ] corrections for each 
energy. The decrease of e(0)5 sv + /hard with energy mainly comes from the efficiency decrease 
(see Table HI). 



VII. SYSTEMATIC UNCERTAINTIES AND RESULTS 



The following have been considered as sources of systematic uncertainty for each con- 
tinuum cross section measurement: luminosity, radiative correction, trigger efficiency for 
hadronic events, multiplicity correction, and hadronic event selection criteria (see Table 

ED). 
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The hadronic vacuum polarization terms in Eq. ([6]) are of the order of 0.01. Since they 
dominate the uncertainty in the radiative correction, we have assigned a 1% uncertainty in 
do from this source. 

Each selection criterion, except for multiplicity, has independently been varied by about 
20% to estimate its systematic uncertainty. Moreover, we used a visible energy E vis /2E heiLm _ 
cut variation of about 40% to assign an additional systematic uncertainty related to the 
two- photon background and the visible energy shift: we assign half the difference between 
R values obtained with cuts on E v i S /2Ei, cam of 0.3 and 0.4 (0.6 and 0.7) as systematic 
uncertainty associated with two-photon background (visible energy shift). 

Even though the multiplicity cut was part of the hadronic event selection criteria, we 
assigned a separate systematic uncertainty to our correction procedure by taking the dif- 
ference between the corrected and uncorrected R values at the nominal cut level. The 
summary of all systematic uncertainties is given in Table IVIl At each energy we divide the 
systematic uncertainty into a common uncertainty that correlated across all energy points 
and an uncorrelated uncertainty that is independent for each energy point. The decrease of 
the uncorrelated systematic uncertainty with decreasing beam energy is mainly due to the 
energy dependence of the two-photon interaction background cross section. 

TABLE VI: Systematic uncertainty (in %) assigned to each energy point. The last two rows show 
common and uncorrelated uncertainties. 



Energy (GeV) 10.538 10.330 9.996 9.432 


8.380 7.380 6.964 


Luminosity 


1.00 


1.10 


1.10 


1.10 


0.90 


0.90 


1.00 


Trigger 


0.09 


0.09 


0.11 


0.08 


0.12 


0.13 


0.19 


Radiative 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


correction 
















Multiplicity 


1.06 


1.38 


0.99 


0.84 


0.43 


0.38 


0.38 


correction 
















Hadron event 


1.51 


1.09 


1.31 


1.31 


1.05 


1.02 


0.79 


selection 
















Total 


2.32 


2.30 


2.21 


2.15 


1.76 


1.74 


1.68 


Common 
Uncorrelated 


1.87 
1.37 


1.67 
1.59 


1.85 
1.22 


1.87 
1.05 


1.62 
0.70 


1.64 
0.57 


1.58 
0.55 



Values of R for each of the energy points are shown in Table IVHI and are presented in 



Fig. [10] along with MARK I [2j, MD1 [22|, Crystal Ball |23j, and prior CLEO [17| results. 
Our new measurements of R are in good agreement with previous MD1, Crystal Ball, and 
CLEO measurements; however, they do not agree with the MARK I results. 

Table IVIIII shows the a s values obtained by solving Eq. ([2]) (with four flavors) for a s 
at each of the seven ener gies (extraction of a s when quark masses are taken into account 
can be found in reference 24]). Comparing a s values with the QCD predictions (25| at our 
energies, which assumes the combined world average of a s (M|) = 0.1189 ± 0.0010, we find 
agreement within our quoted uncertainties. Moreover, we compare our results at 10.330, 
9.996, and 9.432 GeV energies with CLEO measurements @ at T(1S), T(2S), and T(3S) 
energies, which are based on direct photon spectrum studies, and we find a good agreement 
within our quoted uncertainties as well. 



14 



5.5 



5.0 



4.5 



4.0 



3.5 



lit 

iTL 



km 



T □ 



i i 



• CLEO07 

O CLE0 98 

■ Mark I 

A MD1 

□ Crystal Ball 




I.O 8.5 9.0 9.5 

Vs (GeV) 



10.0 10.5 



FIG. 10: R measurement at 5.0 - 10.5 GeV cm. energies. The top figure compares our mea- 
surement with other experiments [20J. The bottom plot shows our new results, where two sets of 
uncertainties represent combined uncorrelated and statistical uncertainties and total uncertainties; 
the line represents Eq. ([2]) with our average A = 0.31 GeV, and the shaded area indicates the 
ii-values corresponding to one standard deviation in the uncorrelated systematic uncertainty in A. 



TABLE VII: Measured values of R(s) with statistical and systematic (common and uncorrelated) 
uncertainties, respectively. 



GeV 


R(s) 


10.538 


3.591 


± 


0.003 


± 


0.067 


± 


0.049 


10.330 


3.491 


± 


0.006 


± 


0.058 


± 


0.055 


9.996 


3.497 


± 


0.004 


± 


0.064 


± 


0.043 


9.432 


3.510 


± 


0.005 


± 


0.066 


± 


0.037 


8.380 


3.576 


± 


0.024 


± 


0.058 


± 


0.025 


7.380 


3.550 


± 


0.019 


± 


0.058 


± 


0.020 


6.964 


3.597 


± 


0.033 




0.057 


± 


0.020 
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TABLE VIII: Derived values of a s (s) with statistical and systematic (common and uncorrelated) 
uncertainties, respectively using the prescription outlined in the text. 



GeV 


a s (s) 


10.538 


0.232 


± 


0.003 


± 


0.061 


± 


0.045 


10.330 


0.142 


± 


0.005 


± 


0.051 


± 


0.049 


9.996 


0.147 


± 


0.004 


± 


0.057 




0.038 


9.432 


0.159 


± 


0.004 


± 


0.058 


± 


0.033 


8.380 


0.218 


± 


0.022 


± 


0.053 




0.023 


7.380 


0.195 


± 


0.017 


± 


0.052 


± 


0.018 


6.964 


0.237 


± 


0.030 


± 


0.052 


± 


0.018 



To test the compatibility with other measurements of a s we note the expected running 
of a s with energy [20| : 



4tt 



a s [s 



/5 ln( S /A 2 ) 



2f3 1 ln[ln(s/A 2 )] 

~M WW) 



+ 



+ 



#ln'( S /A s 
/? 2 /3o 5\ 



x In 



ln(s/A 2 



4 



(14) 



where nj presents the number of quarks which have mass less than y/s/2, A represents 
the QCD energy scale, and the /3-functions are defined as follows: /3 = 11 — 2n//3, fa = 
51 - 19n//3, and /3 2 = 2857 - 5033n//9 + 325^/27. 

To find A, we use our a s values at each energy point and solve Eq. (fl4"]) . assuming rif is 
equal to 4. The value of A varies from 0.11 at 10.330 GeV to 0.67 at 10.538 GeV. Using 
Eq. (fl~4"|) with our average value of A, we extract the value of the a s at yfs = M Z - Our 



results for a s imply A = 0.311^^:21 GeV and a s (Mf ) = 0.126 ± 0.005 where the 

uncertainties represent statistical and total systematic, respectively. 

Our results for a s {M\) and A(n/ = 4) agree with the world averages a s (M z ) = 0.1189 ± 
0.0010 H and A{n f = 4) = 0.29 ± 0.04 GeV @. 

In conclusion, we have shown that our measurements of R(s) are consistent with the 
QCD predictions of Eq. ([1]) and Eq. (TT4"1) with better precision than previous measurements 
of R(s). However, since the Eq. ([!]) dependence of R(s) on a s (s) is only in higher order 
terms, we do not have a determination of a s (s) or A that is competitive with measurements 
by other methods. 
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